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Abstract 

In this paper we present quantitative results for coherent pion production by neutrinos scattered 
off nuclei within the framework developed by Gounaris, Kartavtsev and Paschos. The method is 
based on PCAC and uses helicity cross sections for the scattering of weak gauge bosons on nuclei. 
The process relies on experimental data for elastic pion-nucleus scattering. A detailed analysis of 
the differential and integrated cross sections is presented for neutral and charged currents, with 
special emphasis on the regions of integrations. The results are extended to energies of 10.0 GeV 
and are compared with experimental data. 
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FIG. 1: Feynman diagram of the charged current reaction 
I. INTRODUCTION 

The process of coherent pion production by neutrinos is one of the important reactions 
that occur in neutrino oscillation experiments. It dominates at very low momentum transfers 
of the leptons and in the forward direction which may help to monitor the flux of the beam 
and will help to estimate backgrounds. For instance in the oscillation of into v e one 
searches in the far away detector the reaction v e N —>■ e~X and coherent production of 7r°'s 
produces a background. For this reason there are new theoretical articles and experimental 
mesurements dealing with this process. 

The neutral and charged current reactions are 

v^k)N(jp) - ^(AOAVKV) (1) 
v„(k)N(p) -> ^(k')N(p')n + ( P7T ), (2) 

where N is a nucleus. For the process we use variables in the rest frame of the nucleus with 
q = k — k', Q 2 = —q 2 , v = E — E' and t = (q — p n ) 2 . Coherent production implies that the 
nucleus does not break up or alter its quantum numbers during the process. 

Due to the isospin structure of the charged and neutral currents, the charged cross section 
is at high energies, approximately, twice as big as the neutral one. At low energies differ- 
ence arise from the non- vanishing mass of the outgoing lepton and the Cabibbo-Kabayashi- 
Maskawa-matrix element (CKM) in the charged current case. Hence it is sufficient to con- 
centrate on the charged current reaction within the theory. Its Feynman diagram is shown 
in figure [Q 

Our calculation [1] is based on two basic facts: 
1. Coherent production of pions is the process where the four-momentum-transfer squared 
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between the current and the produced pion is small so the nucleus remains intact. The 
minimum value is given to a good approximation by 

12 , ™2\ 2 



Q 2 + mj 
2u 



(3) 



and is achieved for vM » Q 2 . The important experimental region will turn out to be 
Q 2 < 0.1 GeV 2 . We remind the reader that, neglecting the nucleus recoil energy, t is 
calculated as [2] 
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-I 2 



/i,7T 
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2. In this region, the dominant component of the leptonic current has helicity zero, i.e. 
e^(A = 0) defined in equation below. The domain vM » Q 2 and Q 2 = (a few) m 2 , 
is the region where PCAC is valid and it has been shown that the dominant amplitude is 
determined by chiral symmetry. In fact when we write the amplitude for the subprocess 
WN — > Ntt as the sum 

pion pole + IZ^ 
with 11^ the remainder, the symmetry determines 



%W = -UV2T{tt + N 



7T 



-JV) 



(5) 



which is exactly the term that the zero helicity polarization selects and makes it dominant. 
In this picture the weak current converts into a pion before entering the nucleus. We shall 
quantify these results in the following sections. 



This approach has many advantages: 

1. It avoids discussing what happens within the nucleus since it will use ttN elastic 
scattering data. In this sense it incorporates the pion-nucleus diffrative peak into the 
neutrino scattering. 

2. It is applicable at low and high energies, provided that some required kinematic cuts 
are made. 

3. It includes the lepton mass exactly and quantifies the various approximations. 
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The purpose of this article is to quantify the results with numerical calculations. There 
are other articles on this topic [11 0, Q, E3, E, Is]. One among them follows the same approach 
{si. A very different approach descibes coherent pion production at very low energies as the 
excitation of the A-resonance and includes modifications through nuclear medium effects 



demanding that the nucleus remains in its ground state for coherent scattering 



y,y,QB. 



II. THEORY REVISITED 



The theoretical background of this paper is the same as in For completeness we 
repeat its main features. The invariant matrix element of the charged current reaction is 

GfV U( i 



M 



V2 



-3, (irN\J»\N) 



(6) 



where Gf is the Fermi coupling constant and V u d the CKM matrix element. The leptonic 
current is expressed in a straight forward way 



jp = u(k')% (1 - 75) u(k) 



(7) 



whereas the hadronic matrix element will be treated below. 

The leptonic current is decomposed into the basis of the four polarisation vectors of the 
exchange boson. Select the 3-momentum of along the z-axis and define the basis of 
polarization vectors 
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W 



(8) 



^(A = 0) 



o 



W 



(9) 



Note that they differ from the original paper [1] by their normalization factors. This has 
been taken into account in the following formulas but does not change the results. The 
polarisation vectors satisfy the completeness relation 

•l)V(A)e i/ *(A) 



A=0,±1 



e, e 



I c « 
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(10) 



In the cross sections the leptonic tensor appears 



= 8(kX + - g^k ■ k') + 8ie^k a k ,f3 (11) 

= (12) 

jj' 

with the indices running over the four polarisations. Numerical studies in reference [l| have 
shown that the transverse polarizations, i.e. right and left handed cross sections are smaller 
in coherent scattering and they will be omitted. 

From equation (1121) we obtain the relevant matrix elements: 

\Q 2 (2E„-v) -vmV\ 



2 

m 



Here the mass of the lepton in the final state is explicit and will be kept throughout 
the calculation. We also use later as in reference jl] and the relation between them is 
f —It 



It is now possible to write the spin averaged matrix element as 



G 2 F | V u d | 2 



l-MI 2 = {L^UA? + Lu\ J,ef | 2 + 2L l0 (U)(J^y} . (16) 

As mentioned in the introduction the divergence of the matrix element is determined by 
PCAC. The matrix element of the axial current can be written as a sum of the pion pole 
and the remaining contribution 

~ ^ = ^k T{n+N -> n+N) - (17) 
where T(it + N — > n + N) is the elastic scattering on a nucleus and IZ^ is a smooth function 
which includes all other contributions. The PCAC relation gives 

- W<K = ^rr^ T ^ +N - ^ N )- ( 18 ) 

Combining the two equations one obtains 

<f 7^ = -V2/ 7r r( 7 r + AT -> 7T+A). (19) 



With these relations one is also able to calculate all matrix elements in equation (|T6|) 
because 

a fJl A + 

e l A» + = (20) 

is given by equation f|T8l) and is proportional to ml- 

The matrix element for helicity zero is calculated in the following way. From the property 
e M (A = 0)g^ = it follows that the inner product of the pion pole with the helicity zero 
polarization vanishes. The remaining term is now estimated 

e"(A = 0)7^ = -^=TZ, + O (^pj = ~^J=n- +N - ^ +N ) + O • (21) 

Thus the matrix elements for both polarizations are estimated. 
We write next the triple differential cross section as 



da cc _G 2 F \V ud \ 2 v f 2 J ~ ~ ( m\ \ ~ m 2 1 da n 

^00 + Ml I o +ZL l0^o I o > 17, l^J 



dQ 2 dudt 2(2vr) 2 EIQ 2 \ uu w \Q 2 + m*y Q 2 + m 2 J dt 

where ov is the elastic pion nucleus cross section and where all muon mass terms have been 
kept in the calculation. As mentioned already the right- and left-handed cross sections are 
small and have been neglected. Since they appear as positive additive terms, our estimate 
is a lower bound for the coherent cross section. 

The corresponding neutral current cross section is 

&Qnc G 2 F v fl f da n 

dQ 2 dvdt 4(2tt) 2 El Q 2 00 dt ' 1 ' 

In this formula the muon mass in L Q0 has to be set zero. 



III. METHODS OF INTEGRATION 



The most detailed and convincing evidence for coherent scattering is the explicit obser- 
vation of the triple differential cross sections in equations f l22|) and (T23]) . The characteristic 
signature is the sharp peak in the t-distributions which was the main feature in the original 
discovery and interpretation 0, Q . Subsequent and recent experiments integrate over the 
t- and other variables so that special attention must be given to the ranges of integration 
in order to ascertain that the model is still valid in these regions. In this article we shall 
integrate over data for 7r + C 12 elastic scattering. 
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We identify the incoming pion energy with the variable v and integrate over experimental 



data for elastic pion-nucleus scattering 
follows: 



13] . We parametrize the cross section as 



dt 



aexp[— b\t\] 



(24) 



and fitted the parameters a and b to the data. Their numerical values are given in table [B 



v [GeV] 


a [barn/GeW] 


b [GeV -2 ] 


0.210 


28.526 


159.657 


0.228 


28.659 


147.986 


0.260 


32.012 


129.022 


0.290 


27.162 


101.910 


0.320 


23.600 


90.824 


0.340 


22.734 


90.660 


0.370 


19.000 


83.814 


0.400 


17.924 


84.590 


0.420 


14.594 


73.256 


0.766 


3.759 


49.459 


0.864 


4.172 


58.149 


0.942 


3.649 


56.197 


1.046 


3.523 


53.497 



TABLE I: Parameters of the elastic pion nuclues cross section model. 

As in reference jl| we integrate first the variable t over the range 

/0 2 + m 2 \ 2 , 

< \t\ < oo. 



2u 



(25) 



The upper limit of integration should be the first diflractive minimum and has been expanded 
to infinity because the numerical results are insensitive to values of this cross sections beyond 
the first diflractive minimum, the cross section being already too small. The other limit of 
integration is important at low energies and influences the Q 2 dependence. The integrated 
elastic pion-nucleus cross section 



(26) 
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FIG. 2: Integrated elastic pion nucleus cross section for some momentum transfers. 



depends on Q 2 and v and the results are shown in figure 2. In reference |l| the zv-range was 
limited to 0.4 GeV and now we extended it to v = 1.0 GeV. The data has been extrap- 
olated as constants in regions not covered by data. We observe again a sharp Q 2 -dependence. 



The integration over v should respect the condition 

v » \[Q i (27) 
so that the helicity A = polarization can be expanded to give 

e * (A = 0) " W> + ° (28) 

which has been used in equations (12T1 - 1231) . In order to test this condition we computed 
the double differential cross section d Q% du over the entire range of Q 2 and v using as input 
the curves in figure [2j The results are shown in figures [3] and H] where we also include four 
curves. These four curves correspond to the kinematic regions defined in table [TTJ 

In the earlier article [lJ we used the variable which defined kinematic regions available 
in the experiments. It is a fortunate property of coherent scattering that the cross section 
peaks at low values of Q 2 and v. In fact the contribution to the integrated cross section 
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condition 


value of £ in reference [1] 




t = o 




e = i 




e = 2 


v = 3^02 


e = 3 



TABLE II: Defining regions of phase space. 

from the region 

< v < \[Q i (29) 

is negligible. Our estimates are satisfied for v > l.5y/Q* where most of the cross section 
is located. One can also integrate over restricted regions of the phase space in order to 
determine the fraction of the cross section in these regions. 

00 




FIG. 3: Differential charged current cross section at E v = 1 GeV. The lines represent different 
integration limits (see text). 



The structure of the cross section in figure [3] also indicates the important regions of 
integration. To assure the validity of the approximation in equation (j2Ti) we integrate over 
the range 

max V min \ < V < "max- (30) 

The values for is min and v max are given in the appendix of reference 

The Q 2 integration is treated more carefully. The lower limit of integration is given in 
equation (A8) of reference [1]. In the neutral current case it is zero. 




FIG. 4: Differential charged current cross section at E v = lOGeV. The lines represent different 
integration limits (see text). 

We calculated the differential cross section for E u = 1 GeV and four values for the 
lower limit of the energy v. The results are shown for neutral currents in figure [5] and for 
charged currents in figure [6j The four curves correspond to various cuts in the minimum value 
of the energy u, defined in table[TTJ We note that the cross section between u min < v < \[Q^ 
is negligibly small since the dashed curves for v > v min coincides with curve for v > yQ*. 




FIG. 5: Differential neutral current cross section for £ = 1, 2, 3 (top to bottom) and £ = (dashed) 
at E v = 1 GeV. 
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FIG. 6: Differential charged current cross section for £ = 1, 2, 3 (top to bottom) and £ = (dashed) 
at = 1 GeV. 

We repeated this calculation for various incident neutrino energies and the results are 
shown in figures and M for neutral and charged currents, respectively. The interesting 
feature is that the cross section extends now to higher values of Q 2 . Although the differential 
cross section is concentrated at low momentum transfers, for higher neutrino energies there 
is a tail that extends to larger values of Q 2 and must be taken into account in the integrated 
cross sections. 

The neutral current cross section reaches a specific value at Q 2 = 0, the Adler point. It 
is estimated from our formulas (JTBT) and ( 1231) to give 



da 



NC 



dQ 2 



G 2 F f 2 
2tt 2 



du <t v (v) — 



1 



1 



< ov > < In h — 



2vr 2 



1 



(31) 
(32) 



where < ov > is the weighted elastic pion-nucleus cross section. For E u = 1 GeV the Adler 
point is about 112-10" 40 ^ for < a v >^ 165 mb. This is a realistic average of the weighted 
hadronic cross section (see figure [2]). A similar comparison for the charged current is more 
indirect because of the phase space effects introduced by the mass of the muon. For the 
charged current an integration over a small region of Q 2 is more appropriate for comparison 
with experiments. 
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FIG. 7: Differential neutral current cross section for £ = at E v = 1, 5 and 10 GeV. 
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FIG. 8: Differential charged current cross section for £ = at E v = 1, 5 and 10 GeV. 

It is worth mentioning that for higher energies the turnover of the differential cross section 
in the charged current case is less prominent and disappears, a property caused by the fact 
that the muon mass is negligible compared to the neutrino energy. 

At higher neutrino energies bigger momentum transfers gain more importance. Figures 
[7] and M show the differential cross section of the neutral and charged current process, re- 
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spectively, for £ = 0. For higher energies the cross section is not limited within a small Q 2 
region. Thus the Q 2 integration interval has to be extended to higher values. 



IV. INTEGRATED CROSS SECTION 



The results in figure [8] indicate that the formulas we derived have a tail that extends to 
higher values of Q 2 . However, at large values of Q 2 the nucleus breaks up and the scattering 
is not coherent. At the same time the PCAC approximation does not hold. Fortunately, 
the change in the curves of figures 9 and 20 between £ = 1 and 4 is relatively small. 
Experimental groups separate the coherent events from other events by the absence of stubs 
in bubble chambers, a signature indicating that the nucleus does not break up. The BEBC 
group Ijj] found that the Q 2 -values are lower for coherent events and do not extend beyond 
2.0 GeV 2 . These results suggest that a Q 2 max must be introduced in the integration of the 
cross section or a phenomeno logical factor ^ g 2 "^a J i n order to provide a cut-off in the Q 2 
dependence (see [151]). 



In this work we calculate integrated coherent cross sections for various upper values of 
Qmax = 0-2, 0.5, 1.0 and 4.0 GeV 2 . The results are shown in figure El for neutral current and 
figure [TOl for charged current reactions, where the effect of Q 2 max is important. 
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FIG. 9: Integrated neutral current cross section with Q 
(bottom to top). See 
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max 



0.2, 0.5, 1.0 and 4.0 GeV 2 for £ = 



171 . 1 181 ] for experimental data. 
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FIG. 10: Integrated charged current cross section with Q 



0.2, 0.5, 1.0 and 4.0 GeV 2 for 



£ = (bottom to top). See 
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201 ] for experimental data. 
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or comparison we also included experimental results from several groups 



0, Q, 0, Q, 



201 ] . The experiments use different targets which we must rescale to the carbon target. 



Since the pion-nucleus cross section is a main input in this work we must rescale the results 
according to 



^carbon &exp 



-A-carbon 



(33) 



which has been established in pion-nucleus elastic scattering 



211 ] . With this rescaling the 
agreement at higher energies prefers a Q 2 max between 1.0 and 4.0 GeV 2 . For the energy of 
the neutrinos we use an average value from the neutrino flux. 
There are two ways to account for the Q 2 -dependence: 

1. Whenever possible, introduce an experimental cut-off to the coherent events where 
our model is valid. 

2. Integrate the cross sections on equations (|22|) a nd fl23l) with a phenomenological factor 



introduced by other authors 



10 



15j, which represents the effects of heavier 



vector mesons and treat as a parameter. 
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V. SUMMARY 



Neutrino induced coherent pion production is described adequately by the method de- 
scribed in reference It is argued again that chiral symmetry relates coherent production of 
pions to the pion-nucleus elastic scattering in a general way, provided that Q 2 ~ (a few) -m 2 . 

Our formalism and results are close but not identical to a recent article jsj, where the 
elastic pion-carbon data are now used as well, which reduces their earlier predictions. We 
both use elastic pion-nucleus scattering data and the small difference arise from the handling 
of the experimental data and the limits of integrations. 

Using pion-carbon scattering data we presented quantitative results for many cross sec- 
tions. We emphasize that kinematic limits, as in equation (1251) for t and the maximum value 
of Q 2 in figures ITlfTOl are important. An overview of the double-differential cross section is 
shown in figure [3] and H] showing the main characteristics of the process. The calculation can 
be extended to higher energies when the value of Q 2 max is better understood. 
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